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ABSTRACT 

In an attempt to study whether the formation of brown dwarfs (BDs) takes place as a 
scaled-down version of low-mass stars, we conducted IRAM30m/MAMBO-II obser- 
vations at 1.2 mm in a sample of 12 proto-BD candidates selected from Spitzer/IRAC 
data in the B213-L1495 clouds in Taurus. Subsequent observations with the CSO at 
350 /mi, VLA at 3.6 and 6 cm, and IRAM 30 m/EMIR in the 12 CO(1-0), 13 CO(l- 
0), and N2H" 1 " (1-0) transitions were carried out toward the two most promising 
Spitzer/IRAC source(s), J042118 and J041757. J042118 is associated with a compact 
(< 10 arcsec or < 1400 AU) and faint source at 350 //m, while J041757 is associated 
with a partially resolved (~ 16 arcsec or ~ 2000 AU) and stronger source emitting 
at centimetre wavelengths with a flat spectral index. The corresponding masses of 
the dust condensations are ~ 1 and ~ 5 Mj up for J042118 and J041757, respectively. 
In addition, about 40 arcsec to the northeast of J041757 we detect a strong and 
extended submillimetre source, J041757-NE, which is not associated with NIR/FIR 
emission down to our detection limits, but is clearly detected in 13 CO and N2H + at 
~ 7 km s -1 , and for which we estimated a total mass of ~ 100 Mj up , close to the 
mass required to be gravitationally bound. In summary, our observational strategy 
has allowed us to find in B213-L1495 two proto-BD candidates and one pre-substellar 
core candidate, whose properties seem to be consistent with a scaled-down version of 
low-mass stars. 

Key words: Stars: brown dwarfs - Stars: formation - Stars: individual objects: 
SSTB213 J041757, SSTB213 J042118 - ISM: molecules - ISM: dust - ISM: clouds 



1 INTRODUCTION 
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The formation mechanism of brown dwarfs (BDs) is cur- 
rently a matter of vigorous debate. Since the typical thermal 
Jeans mass in molecular cloud cores is around 1 Mq , a cloud 
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should not fragment in cores of substellar masses and the 
formation of brown dwarfs cannot be directly explained as a 
scaled-down version of low-mass star formation. A possible 
scenario is that BDs form from fragmentation of turbulent 
cores. In this scenario, turbulence produces particular den- 
sity enhancements, decreasing the Jeans mass and allowing 
the formation of cloud cores of very low masses (Padoan & 
Nordlund 2004). Another possible scenario is the formation 
of BDs through fragmentation of massive discs and subse- 
quent ejection (Rice et al. 2003; Stamatellos & Whitworth 
2009) , or ejection could also take place in the very first stages 
of formation of a multiple system (e.g., Reipurth & Clarke 
2001; Umbreit et al. 2005; Basu & Vorobyov 2012), depriving 
the ejected objects from gaining further mass. Other scenar- 
ios propose that an initially typical core is photoerosioned 
due to the nearby presence of a massive star (Whitworth & 
Zinnecker 2004). 

An observational way to distinguish among these dif- 
ferent scenarios is the search and characterization of BDs 
in their very first evolutionary stages, what we call here 
the proto-BD stage (following the nomenclature of Pound 
& Blitz 1993; 1995, and equivalent to Class 0/1 stages of 
low-mass star formation). If BDs form in a similar way 
as their low-mass stellar counterparts, one would expect 
to find proto-BDs surrounded by substantial discs and en- 
velopes, with the envelopes similar to the ones observed in 
the first stages of low-mass protostars (e.g., Andre, Ward- 
Thompson, & Barsony 1993), which are strongly emitting 
in the millimetre/submillimetre range. In addition, proto- 
BDs should be found associated with outflow and ejection 
phenomena, whose parameters should be a scaled-down ver- 
sion of low-mass stars, and perhaps should even form plan- 
etary systems. As a matter of fact, a significant number of 
BDs with ages of a few Myr have been shown to have discs 
(e.g., Luhman et al. 1997; White et al. 1999; Fernandez & 
Comeron 2001; Barrado y Navascues & Martm 2003; Pas- 
cucci et al. 2003; Barrado y Navascues et al. 2004, 2007; 
Natta et al. 2004; Apai et al. 2005; Mohanty et al. 2005; 
Jayawardhana et al. 2006a,b; Luhman et al. 2005, 2007, 
2008, 2010; Harvey et al. 2010; Monin et al. 2010; Joergens 
et al. 2012; Riaz et al. 2012b), but most of these discs were 
inferred from optical or near/mid-infrared data, being thus 
T-Tauri analogs or Class II sources. 

Different research groups have searched for proto-BDs 
in the millimetre/submillimetre range. One of the pioneer- 
ing works was that of Pound & Blitz (1993; 1995), who find 
no clear-cut evidences of proto-BDs. These works were fol- 
lowed up by Greaves et al. (2003), who survey an area of 
the Ophiucus star-forming region with the JCMT at 850 /im, 
and find a dozen of possible planetary-mass isolated objects. 
However, their submillimetre sources could be tracing either 
truly proto-BDs or transient objects, and these objects re- 
main to be confirmed through molecular line observations. 
Other candidates to be BDs in the Class 0/1 phase are the 
so-called Very Low Luminosity Objects (VeLLOs), objects 
embedded in dense cores with internal luminosities < 0.1 Lq 
(Di Francesco et al. 2007; Dunham et al. 2008), such as 
L1014-IRS and L1148-IRS, which are associated with dusty 
and dense envelopes, centimetre emission (L1014: Shirley 
et al. 2007), and CO (1-0) and CO (2-1) outflows (Bourke 
et al. 2005; Kauffmann et al. 2011), all properties sugges- 
tive of accretion/ejection of matter. However, for L1014-IRS 



and the other few VeLLOs studied in detail (e.g., Andre, 
Motte, & Bacmann 1999; Young et al. 2004; Kauffmann et 
al. 2005, 2008, 2011; Lee et al. 2009; Dunham et al. 2010) 
their true substellar nature is not clear (eg., Lee 2007), nor 
from the bolometric luminosity (> 0.05 Lq) neither from 
the mass of the envelope, which could increase the present 
mass well above the substellar limit once accreted on to 
the object (Kauffmann et al. 2011). On the other hand, the 
studies of Klein et al. (2003), and Scholz et al. (2006), who 
carry out surveys with the JCMT and/or the IRAM30m 
Telescope, were addressed towards samples of relatively old 
brown dwarfs (> 1 Myr), selected from optical/infrared sur- 
veys, while the most embedded objects are expected at ages 
< 0.5 Myr. Finally, individual proto-BD candidates studied 
in detail in the millimetre/submillimetre range either are 
more evolved than Class I sources (e.g., Phan-Bao et al. 
2008; Harvey et al. 2012; Riaz et al. 2012a), or do not have 
the mass well constrained and could still be stellar (e.g., 
Scholz et al. 2008). 

In this paper we show the results of a millime- 
tre/submillimetre study of a sample of 12 proto-BD candi- 
dates selected in the B213-L1495 dark clouds of the Taurus- 
Auriga complex using Spitzer data. We show the results of 
the observations at 1.2 mm toward this sample, as well as 
submillimetre data for the two most promising candidates 
found at 1.2 mm. One of the candidates was not previously 
reported, while the other candidate was already reported in 
Barrado et al. (2009), and here we present subsequent ob- 
servations aimed at further constraining its nature, which 
additionally reveal a nearby starless core of substellar mass, 
i.e., a possible pre-substellar core. 



2 SELECTION STRATEGY AND 
OBSERVATIONS 

2.1 Selection strategy 

We used IRAC and MIPS (Spitzer) data toward one of the 
youngest regions in Taurus, the B213-L1495 region. This 
region was selected from Froebrich (2005) list of Class 
sources, because it was one of the nearby (< 300 pc) regions 
with a higher ratio of Class 0/1 to Class II IRAC sources 
(see Fig. 1 in Barrado et al. 2009). For B213-L1495, we 
downloaded the Spitzer archive data and selected 12 IRAC 
sources which were classified as Class 0/1 according to IRAC 
color-color diagrams (Allen et al. 2004), and substellar (with 
a magnitude below the 0.075 Mq cut in the color-magnitude 
diagram, after comparing with evolutionary models, Bar- 
rado et al. 2009). In addition, we rejected those objects 
classified as extragalactic candidates following the criteria 
of Harvey et al. (2006), Joergensen et al. (2006), and Bouy 
et al. (2009). Fig. Q] shows the final sample overlaid on an 
extinction map of the region, elaborated for this work using 
the star count method on 2MASS data. 
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Figure 1. Extinction map in 2MASS J-band derived with 
an adaptive star count method (Cambresy et al. 1997) in the 
B213-L1495 region (45' radius). Members and candidate mem- 
bers from Kenyon &; Hartmann (1995) and references therein 
have been taken into account (removed from the star count- 
ing) when computing the Aj values. The resolution of the map 
varies from 30—50 arcsec, depending on the density of stars de- 
tected by 2MASS. Plus signs mark the positions of the selected 
Spitzer/IRAC sources observed with the IRAM30m Telescope, 
and the stars correspond to the sources additionally observed with 
the CSO. 



2.2 IRAM 30 m ON-OFF at 1.2 mm and APEX at 
870 (ira for the sample of 12 proto-BD 
candidates 

We observed with the IRAM 30 rr0 single-dish Telescope 
(Pico Veleta, Spain) the 1.2 mm continuum emission of the 
12 selected BD candidates using the 117-channel MAMBO- 
II bolometer array. The Half Power Beam Width (HPBW) of 
the telescope at 1.2 mm is 11 arcsec. Observations were car- 
ried out in 2007 May, October, November, and 2008 Novem- 
ber 14, and the opacities ranged between 0.30 and 0.40. Each 
source was centred in the most sensitive pixel of the array 
(pixel 20). Typical pointing corrections were < 5 arcsec, pro- 
viding a positional accuracy of ~ 3 arcsec. The background 
subtraction was performed by using the ON-OFF observ- 
ing mode, with a wobbler throw of 32 arcsec. Every source 
was observed for at least one scan of 20 minutes with an 
integration time per subscan of 60 s. The data were reduced 
using the MOPSIC pipeline provided by IRAM. In order to 
test the consistency of the data we measured for each source 
the flux density for each scan separately, and obtained val- 
ues that were consistent within the uncertainties. In Table [T] 

1 The IRAM 30m Telescope is operated by the Institut de Ra- 
dioastronomie Millimetrique, which is supported by INSU/CNRS 
(France), MPG (Germany), and IGN (Spain). 



we summarize the observations presented in this work, and 
in Table [2] we give the list of the 12 selected Spitzer/IRAC 
sourcefl total integration time, opacity at 250 GHz, and 
the measured flux density at 1.2 mm. The final rms noise at 
1.2 mm was typically 0.7-1.5 mjy beam -1 , depending on the 
integration time. None of the 12 Spitzer/IRAC sources was 
detected above the 4cr level, corresponding to envelope/disk 
masses < 3 Mj up (for a dust temperature of 15 K and 
a dust opacity at 1.2 mm of 0.009 cm 2 g _1 , Ossenkopf & 
Henning 1994, Table [2} • However, for two of the observed 
Spitzer/IRAC sources we measured fluxes above 2a, namely 
J041757 and J042118, which were subsequently observed 
with the CSO. 

In addition, we used public archived data from 
LABOCA for the Atacama Pathfinder Experiment (APEX) 
antennc0. Observations at 870 fim were carried out in 
2009 July under the project E-083.C-0453A.2009 (P.I. M. 
Tafalla). The region mapped corresponds to Taurus B213 
and B211 filaments, and was observed during a total time of 
14 hours, under good weather conditions (precipitable water 
vapor (pwv) ranged between 0.25 and 1.8 mm). Data reduc- 
tion was performed using BoA and miniCRUSH software 
packages following standard procedures. The LABOCA im- 
age covers 5 of the 12 Spitzer/IRAC sources of our sample, 
and none of the 5 sources were detected (see Table [2] for 
upper limits). 



2.3 CSO at 350 fim. for J042118 and J041757 

The Caltech Submillimetre Observatory (CSO) observations 
toward J042118 and J041757 were carried out on 2008 Jan- 
uary 28, using the SHARC-II bolometer array at 350 fim. 
The HPBW of the CSO at this wavelength is 8.6 arcsec, 
which was convolved to a final beam of 10.4 arcsec. Typi- 
cal opacities at 225 GHz were around 0.04. The observations 
were carried out in imaging mode using a lissajous pattern of 
20 arcsec of amplitude and 14.142 seconds of period, result- 
ing in a final image size of 3.3 x 3.2 arcmin 2 . In this observing 
mode, the telescope modulates the X and Y position with a 
different sine wave for each axis (starting from the centre of 
the map). As a consequence, the rms is lower in the centre 
of the image than in the edges. Each scan was integrated 
for 10 minutes, and the total observing time per source was 
around 1-1.5 hours. The absolute flux uncertainty was esti- 
mated from a measurement of the flux of HLTau, and was 
< 10%. The uncertainty in position was ~ 3 arcsec. The 
reduction of the SHARC-II data was carried out with the 
CRUSH software (Kovacs 2008), version 1.62-4. For both 
sources we assumed that the emission is compact (< 30 arc- 



2 Throughout this paper, we will refer to each source of Table [2] 
designated by 'SSTB213 J04mmss.ss+ddmmss.s' as 'J04mmss'. 

3 APEX is a collaboration between the Max-Planck-Institut fuer 
Radioastronomie, the European Southern Observatory, and the 
Onsala Space Observatory. 
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Table 1 . Summary of observations and catalogs used in this work 





Wavelength 


Observing 


Sources 


Date 


HPBW a 


P. A. a 


Telescope/Instrument 


(jjzd) 


mode 


observed 


(yy/mm/dd) 


(arcscc) 


(°) 


CFHT / McgaPrimc-McgaCam 


0.75,0.90 


photometry 


J041757 


04/12/05 b 


~ 1 




APO/SDSS 


0.89 


photometry 


J042118 


> 2004 


~ 1 




CAHA/Omega2000 


1.25,1.65,2.17 


photometry 


J041757 


07/ll/04 b 


~ 1 




UKIDSS/WFCAM 


1.25,1.63,2.20 


photometry 


J042118 


07/09/28 


~ 1 


— 


Spitzer/IRAC 


3.6,4.5,5.8,8.0 


photometry 


Table 2 


05/02/21 b 


1.4-1.9 




WISE/SC 


12,22 


photometry 


J041757,J042118 


2010 


6.5-12 




Spitzer/MIPS 


24,70 


photometry 


Table 2 


05/02/28 b 


6-18 




CSO/SHARC-II 


350 


cont. imaging 


J041757,J042118 


08/01/28 


10 




APEX/LABOCA 


870 


cont. imaging 


5 sources Table 2 


09/07 


18 




IRAM 30m/MAMBO-II 


1200 


ON-OFF 


Table 2 


07/05,07/11 


11 




IRAM 30m/EMIR 


2700, 13 CO(1-0) 


line imaging 


J041757 


09/10/28 


22 




IRAM 30m/EMIR 


2600, 12 CO(1-0) 


line imaging 


J041757 


09/10/28 


22 




IRAM 30m/EMIR 


3200,N 2 H+(l-0) 


freq. switch. 


J041757 


10/08/13 


27 




VLA 


36000 


cont. imag.,D c 


J041757 


08/09/04 


10 x 9 


-68° 


VLA 


36000 


cont. imag.,B c 


J041757 


09/02/25 


1.7 x 1.6 


-89° 


VLA 


60000 


cont. imag.,D c 


J041757 


08/09/04 


17 x 16 


+77° 


VLA 


60000 


cont. imag.,B c 


J041757 


09/02/25 


2.4 x 1.8 


+61° 



Half Power Beam Width (HPBW) and Position Angle (P. A.) of the beam. 
Some data already published in Barrado et al. (2009). 
'B' or 'D' stands for the array configuration. 

Table 2. IRAM 30m/MAMBO-II at 1.2 mm and Archive APEX/LABOCA at 870 fira observations for the sample of proto-BD candidates 



Position 



Source 



q(J2000) 



<5(J2000) 



Int. time a 
(min) 



. 2mm b 

(mjy) 



A/f c r>870fjm d 
(M Jup ) (mjy) 



SSTB213 J041726.38+273920 


.0 


04 


17 


26 


.38 


+27 


39 


20 


.0 


40 


< 


4.2 


< 


2 


.7 


< 


84 


SSTB213 J041740.32+282415. 


.5 


04 


17 


40 


.32 


+28 


21 


15 


.5 


40 


< 


4.2° 


< 


2 


.7 






SSTB213 J041757.77+274105. 


.0 


04 


17 


57 


.77 


+27 


41 


05 


.0 


100 


< 


2.9 b 


< 


1. 


.8 


< 


84 


SSTB213 J041828. 08+274910. 


.9 


04 


18 


28. 


.08 


+27 


49 


10 


.9 


20 


< 


5.4 


< 


3 


.4 


< 480 


SSTB213 J041836. 33+271442. 


.2 


04 


18 


36. 


.33 


+27 


14 


42. 


.2 


20 


< 


7.2 


< 


4. 


.6 


< 


96 


SSTB213 J041847.84+274055. 


.3 


01 


18 


47. 


.84 


+27 


40 


55. 


.3 


20 


< 


5.8 


< 


3 


.7 


< 


88 


SSTB213 J041913. 10+274726. 


.0 


01 


19 


13. 


.10 


+27 


47 


26. 


.0 


20 


< 


5.9 


< 


3 


.7 






SSTB213 J041938.77+282340 


.7 


04 


19 


38 


.77 


+28 


23 


40 


.7 


20 


< 


5.0 


< 


3 


.2 






SSTB213 J042016.70+280033 


.7 


04 


20 


16 


.70 


+28 


00 


33 


.7 


20 


< 


5.7 


< 


3 


.6 






SSTB213 J042019.20+280610 


.3 


04 


20 


19 


.20 


+28 


06 


10 


.3 


20 


< 


5.4 


< 


3 


.1 






SSTB213 J0421 18.43+280640 


.8 


04 


21 


18 


.43 


+28 


06 


10 


.8 


10 


< 


4.0 b 


< 


2 


.5 






SSTB213 J042123.70+281800 


.6 


04 


21 


23 


.70 


+28 


18 


00 


.6 


20 


< 


5.6 


< 


3 


.5 







Total integration time for IRAM 30m/MAMBO-II ON-OFF observations. 

Upper limits correspond to 4o\ For the case of J041757 and J042118, the weighted mean of the different subscans was 1.89 ± 0.72 and 
2.15+ 1.00 mjy respectively, both above 2a, while for the other sources the weighted mean was well below 2a. 

Masses of disc/envelope estimated from the 1.2 mm continuum emission, assuming the emission is optically thin, a dust temperature of 
15 K, and a dust opacity at 1.2 mm of 0.009 cm 2 g" 1 (Ossenkopf & Henning 1994). 

Upper limits correspond to 4a, and are measured in the Archive APEX/LABOCA image of the B213/B211 filament (see ij 12. 21 1. 

No reliable flux due to the presence of a submillimetre source (J0417426+282442, detected by AKARI with a flux of 3.6 Jy at 140 fim) 
within the wobbler throw. We estimated that up to 6 mjy could have been subtracted to this source. This is the only source of our 
sample with an AKARI submillimetre source within the wobbler throw of the IRAM 30m/MAMBO observations. 



sec and faint (< 1 Jyfl. The rms noises achieved at the 
map centres are 9 mjy for J042118 and 6 mjy for J041757. 

4 In Barrado et al. (2009) we show the CSO emission for J041757 
assuming that the source is extended, i.e., allowing to recover 
emission from structures > 30 arcsec. For this case the source 
associated with J041757 is ~ 30 arcsec, and has a small exten- 



sion to the north-east. However, the 'compact' option used in this 
work, although yielding more negative features, recovers a struc- 
ture which is more consistent with independent observations of 
the region (as shown in this paper). For the J042118 region, the 
assumption of 'extended' emission (forcing the software to recover 



2.4 VLA observations at 6 and 3.6 cm for 
J041757 

We obtained radio continuum data at 6 cm (4.86 GHz) 
and 3.6 cm (8.44 GHz) towards the source J041757 us- 
ing the Very Large Array (VLA) of the National Radio 
Astronomy Observatory (NRAOjj. Observations were per- 
formed on 2008 September 4 (project AP566), and on 
2009 February 25 (project AD593), with the array in D 
and B configuration respectively. We used the standard 
VLA continuum mode (4IF, 50 MHz per IF), and both 
right and left polarizations were processed. The phase cen- 
tre of the observations was R.A.(J2000.0)=04 h 17 m 57.8 3 , 
Dec(J2000.0)=+27°4l'05.0" in both epochs. For the first 
set of data (D configuration of the array) and for each fre- 
quency, we used ~ 40 minutes on source, and for the second 
set (B configuration) we integrated ~ 1.5 hour on source at 
each frequency. See Table [3] for further details on the obser- 
vations. 

Both sets of data were obtained during the Extended 
Very Large Array (EVLA) transition mode. Calibration and 
data reduction were performed using the Astronomical Im- 
age processing System (AIPS, developed by NRAO), and 
following NRAO official recommendations for the EVLA 
transition mod«[f|. As a result, data in D configuration at 
both 6 and 3.6 cm were successfully calibrated. As for the 
second dataset (B configuration), anomalous data affecting 
the flux calibrator were flagged at 6 Thus, we used 

data in B configuration just to illustrate the morphology 
of the emission at higher angular resolution, and for flux 
measurement purposes we used only the first dataset (D 
configuration). For the D configuration data, we produced 
radio continuum maps with natural weighting to optimize 
the sensitivity of the images, achieving rms noises of ~0.025 
and ~0.055 mjy beam -1 at 6 cm and 3.6 cm, respectively. 
For the B configuration data, we produced maps with nat- 
ural weighting and using a taper in the uv-d&ta, of 150 kA 
at 6 cm and 120 kA at 3.6 cm. For the B configuration data 
we refrain from giving the rms noises because of anomalous 
data for the flux calibrator. 

structures < 60 arcsec to avoid sky gradients) yields very similar 
results to the assumption of 'compact' emission. 

5 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 

6 During the transition to the EVLA, the 50 MHz continuum 
mode was affected by high closure errors on EVLA- VLA baselines 
(affecting strongly the 3.6 cm receiver; for this reason, Barrado 
et al. (2009) report only the flux density at 6 cm). This problem 
is caused by non-matched bandpass shapes of different type of 
antennas. In this work, we tried to fix this problem by observing 
a strong source with known structure and determining baseline- 
based closure corrections using the AIPS task BLCAL. 

7 For the B configuration data, anomalous jumps in the ampli- 
tudes and phases for certain EVLA antennas, specially during 
the end of the target and flux calibrator observation, prevented 
us from a reliable flux calibration. These anomalous data were 
flagged only for the 6 cm case, because the 3.6 cm observations 
were too short. However, flagging these data with anomalous 
jumps improved the signal-to-noise, while the morphology of the 
object was found unaffected, giving us confidence on the morphol- 
ogy of the source seen in the B configuration data at both 6 and 
3.6 cm. 



A search for proto-BDs in Taurus 5 

Table 3. Parameters of the VLA radio continuum observations 
towards J041757 



A Flux Flux Phase Flux a 
Config. (cm) calib. (Jy) calib. (Jy) 



D 6 3C147 7.9 J0431+206 2.34 ± 0.03 

D 3.6 3C147 4.7 J0431+206 1.590 ± 0.005 

B 6 3C138 3.6 J0431+206 2.63 ± 0.05 

B 3.6 3C138 2.4 J0431+206 1.06 ± 0.04 

a Bootstrapped flux. 



2.5 IRAM 30m CO, "CO and N 2 H+ emission for 
J041757 

J041757 was observed in the CO (1-0) (115.271204 GHz), 
and 13 CO(1-0) (110.201370 GHz) transitions with the 
IRAM 30 m Telescope using the EMIR (Carter et al. 2012) 
heterodyne receivers. The observations were carried out on 
2009 October 28 and 29, with system temperatures of 200- 
400 K (r ~ 0.16-0.35). We conducted on-the-fly maps in 
position-switching mode, scanning the region in R. A., and 
with steps of 4 arcsec in Dec, covering a size of 1.6 x 

1.6 arcmin 2 centred on J041757. For each of the CO (1-0), 
and 13 CO(1-0) transitions we used two units of the VESPA 
(VErsatile SPectrometer Array) backend, one in each polar- 
ization, of 40 kHz (0.11 km s -1 ) of spectral resolution and 
a bandwidth of 40 MHz (~ 100 km s _1 ). Simultaneously to 
the VESPA backend, we also used the WILMA (Wideband 
Line Multiple Autocorrelator) backend, providing a band- 
width of 9000 km s _1 with ~ 5 km s" 1 of spectral resolution 
at the frequency of the CO (1-0) transition. The OFF posi- 
tion used, at (—7600, —287) arcsec from the map centre, was 
selected as free of CO from a survey in Taurus (Goldsmith 
et al. 2008), and was proved to be free of emission. Cali- 
bration following the standard procedures was carried out 
during observing time and using the MIRA package, and 
typical calibration errors are within 15-20%. To make the 
final map, spectra were averaged every 8 arcsec in both RA 
and DEC. Typical rms noise of spectra of each point of the 
final map were around 0.24, and 0.12 K, for CO (1-0), and 
13 CO (1-0), respectively, in main beam brightness temper- 
ature scale. The main beam efficiency at the CO (1-0) line 
frequency is 0.82 (F cff =94.01%, B cff =77.53%) and at the 
13 CO(1-0) frequency is 0.83 (F cff =94.18%, B cff =78.13%). 

Subsequent observations of the dense gas tracer 
N 2 H+(l-0) (at 93.173764 GHz) were carried out on 2010 Au- 
gust 13 under reasonable summer weather conditions (pwv 
~ 4 mm; system temperature ~ 130 K; opacity ~ 0.4), in 
frequency switching mode, and using VESPA with 40 MHz 
of bandwidth and 0.02 MHz of spectral resolution. We inte- 
grated at different positions in and around J041757, sep- 
arated ~ 20 arcsec and covering offsets with respect to 
J041757 from (40,40) arcsec to (-20,-20) arcsec, reach- 
ing rms noises in each individual spectra of ~ 0.015 K 
(main beam brightness temperature). The main beam ef- 
ficiency at the N 2 H+(l-0) frequency is 0.85 (F eS =94.76%, 
B c ff =80.15%). All data reduction was done using the 
GILDAS software package supported at IRAM. 
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3 RESULTS 



3.1 J042118 



Although the 1.2 mm emission for J042118 was very 
marginal (2.2<r, Table [2}, the CSO observations at 350 fim 
reveal one clear source with a flux level of 5a ~ 46 mjy 
(Fig. [2]-top). The position of the submillimetre source, 
which is unresolved, is coincident with the position of the 
Spitzer/IRAC and UKIDSS source, and there is marginal 
extended and faint emission at 350 /im about 10 arcsec to 
the west. In Table [4] we give the position, size, flux den- 
sity and estimated mass of the submillimetre source de- 
tected in J042118. The mass of the disc/envelope obtained 
for J042118 is in the range 0.3-3 Mj up , assuming a distance 
of 140 pc, a dust temperature in the range 10-20 K, and 
a dust opacity law at 350 /im of 0.101 cm 2 g _1 (Ossenkopf 
& Henning 1994). Such a mass is at least 3 orders of mag- 
nitude larger than the disc masses estimated from Herschel 
observations for a small sample of young BDs (Harvey et al. 
2012). 



3.2 J041757 

3.2.1 CSO 350 (im and APEX 870 fj,m continuum 
emission 

J041757, with a 2.6a flux at 1.2 mm (Table 0), is well de- 
tected at 350 /mi, and presents a complex emission, with 
one partially resolved source (~ 16 arcsec of size) associ- 
ated with the selected IRAC source J041757, at the centre 
of the field of view, and one stronger and more extended 
source (~ 20 arcsec) to the north-east (hereafter, J041757- 
NE), which is elongated in the southwest-northeast direction 
(Fig. [2}bottom and Table [4]). The flux density of J041757 
at 350 /im is 165 mjy, yielding a mass of the envelope of 
1-10 Mj up , (for a dust temperature of 10-20 K). Assum- 
ing a mean molecular weight (of H2 molecule) of 2.8 mH, 
and a size (diameter) of 2200 AU, we obtain a density of 

< 2.1 x 10 s cm" 3 , a H 2 column density of < 4.6 x 10 21 cm" 2 , 
and A v < 4.9 mag (following Frerking et al. 1982). Sim- 
ilarly, for J041757-NE we estimate a mass of 2-20 Mj up , 
a density of < 2.4 x 10 5 cm" 3 , a H2 column density of 

< 6.4 x 10 21 cm" 2 , and A v < 6.7 mag. These A v val- 
ues are in perfect agreement with the values derived by 
Schmalzl et al. (2010) for this region. It is interesting to 
note that J041757-NE is also recovered in the archived 
APEX/LABOCA 870 /mi image (see Fig. rjjbottom), for 
which we measured a peak intensity of 100 mjy beam" , 
and a flux density (within the 2a contour) of 420 mjy. This 
corresponds to a mass of 75-160 Mj up (assuming a dust tem- 
perature of 10-15 K, and a dust opacity law at 870 /im of 
0.018 cm 2 g-\ Ossenkopf & Henning 1994). The mass of the 
J041757-NE source as estimated from the APEX/LABOCA 
image is about a factor of 8 higher than the mass estimated 
from the CSO because structures > 30 arcsec were filtered 
out during the reduction process of the CSO data (§ 12. 3p . 
Thus, the total mass of J041757-NE is ~ 100 M Jup , from 
which ~ 10 Mj U p are found in a compact structure. 
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Figure 2. CSO 350 /tm signal-to-noise contours superimposed 
on a NIR image. Top: J042118. Contours are —4, —3 (dashed), 
3, 4, and 5 (solid) times the rms noise of the map (noise is 
variable across the map, from 9 (centre) to 11 mjy beam -1 
(wedges)). Greyscale is the UKIDSS i^-band image. Bottom: 
J041757. Black contours correspond to the CSO 350 /im emission 
and are -5, -4, -3 (dashed), 3, 4, 5, 6, 7, 8, 9 and 10 (solid) 
times the rms noise of the map (noise is variable across the map, 
from 6 (centre) to 9 mjy beam" (wedges)). White contours cor- 
respond to the APEX 870 /im emission and are 1, 2, 3, 4, and 5 
times 19 mjy beam . The APEX beam at 870 /im is 18 arcsec. 
Greyscale is CAHA/Omega2000 H-band image (see Barrado et 
al. 2009 for details on the observations), showing the two NIR 
sources, A and B, found near the IRAC source. In both panels, 
the CSO beam, shown in the bottom-right corner, is 10 arcsec, 
and the crosses indicate the position of the Spitzer/IRAC source. 



3.2.2 VLA 3.6 and 6 cm continuum emission 

The VLA centimetre emission around J041757 reveals one 
single source at both 3.6 and 6 cm wavelengths at the po- 
sition of J041757, with no significant emission at the posi- 
tion of J041757-NE. Using D-configuration only, the emis- 
sion is unresolved at both frequencies, and it is detected at 
~7<r at 6 cm, and marginally detected at ~4<r at 3.6 cm, 
with a flux density at both wavelengths of ~ 0.2 mjy (Ta- 
ble [5]). The spectral index in the range of 6 to 3.6 cm is 
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Table 4. Parameters of the sources detected with CSO/SHARC-II at 350 





rms a 


Position 10 


Deconvolved 


Deconvolved 


Deconvolved 


rpeak c 

1 V 


si 


Mass d 


Source 


(mjy) 


o(J2000) 


<5(J2000) 


ang.size (arcsec) 


P. A. (°) 


size (AU) 


(mjy beam -1 ) 


(mjy) 


(M Jup ) 


J042118 


8.9 


04:21:18.38 


28:06:40.1 


10.5 X 8.7 


58.8 


1500 x 1200 


46 ± 9 


16 


0.3-3 


J041757 


6.1 


04:17:57.67 


27:41:05.8 


19.5 x 12.7 


-67.3 


2700 x 1800 


37 ±6 


165 


1-10 


J041757-NE 


8.6 


04:18:00.30 


27:41:36.3 


29.2 x 15.4 


29.7 


4100 x 2100 


84 ±9 


356 


2-21 



a Rms noise at the position of the source. Total integration time was around lh 30min for J041757 and around lh for J042118. 
b Positions corresponding to the intensity peak. 

c Peak intensity and flux density derived by computing statistics within the aperture corresponding to the 3(7 level of the source. 
d Masses derived assuming a dust temperature of 10-20 K, and a dust mass opacity coefficient from Osscnkopf & Hcnning (1994, see 
main text). The uncertainty in the masses due to the opacity law is estimated to be a factor of 2. 



Table 5. VLA radio continuum emission detected at 3.6 and 6 
cm towards J041757 



A Position a Size, P.A. a 

(cm) a(J2000) <5(J2000) (mjy) (arcsec 2 , °) 



6 04:17:57.74 +27:41:04.2 0.20±0.05 2.4 X 1.3, 45 
3.6 04:17:57.75 +27:41:04.4 0.19+0.11 2.9 X 0.8, 40 

a Results based on observations in B configuration of the VLA 
(see § 12. 41 1. and derived from Gaussian fits. Uncertainties in sizes 
are around 0.5 arcsec and in position angles (P. A.) are around 
10°. Sizes and P. A. are deconvolved. 

b Results based on observations in D configuration of the VLA. 
Uncertainties are 2c. 

06-3.6 cm=0.0±0.5 {la error; where S v oc u a ), which is in- 
dicative of optically thin free-free emission (typically aris- 
ing in shocks for the case of low-mass young stellar ob- 
jects and VeLLOs, see § I5.2.3[) . To extract the morpho- 
logical information of the centimetre emission, we imaged 
the B-configuration data only, and the resulting images are 
presented in Figure [3] The figure shows first that the cen- 
timetre source is clearly associated with the NIR source 'B' 
detected by Barrado et al. (2009), while source 'A' shows 
no centimetre emission. Furthermore, the figure reveals at 
both 6 and 3.6 cm a faint extension in the southwest- 
northeast direction. A Gaussian fit to the source yields, 
at both wavelengths, a deconvolved size of approximately 
~ 2.5 x 1.0 arcsec 2 , oriented at PA~ 40 ° (Table 0. 

3.2.3 IRAM 30m molecular emission 

The CO (1-0) spectra averaged over a region of 40 arcsec of 
diameter around J041757 (left) and J041757-NE (right) are 
shown in Fig.[4]-top. The spectra reveal the three well-known 
velocity components associated with the Taurus molecular 
cloud, at around 5, 7 and 8 km s _1 (e. g., Heyer et al. 1987; 
Mizuno et al. 1995). Among the three main velocity compo- 
nents seen in Taurus, the velocity maps of Fig. 6 of Mizuno 
et al. (1995) show that at the position of J041757 there are 
two intersecting filamentary structures with different veloc- 
ities: one in the 4.5-6.0 km s _1 range, and the other in the 
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Figure 3. High angular resolution VLA radio continuum emis- 
sion (red contours) towards J041757, obtained with the array in 
B configuration, superimposed on the Omega2000 //-band image 
(greyscale) and the CSO 350 /im image (black contours). Top: 
6 cm emission. Red contours are —4, —3 (dashed), 3, 4, 6, 8, 
10, and 12 (solid) times the rms of the map (rms not given due 
to technical problems with the flux density scale during the B- 
configuration observations, see § 12.411 . The synthesized beam is 
shown in the bottom-right corner, and is 2.45 X 1.77 arcsec 2 , with 
PA=61°. Bottom: 3.6 cm emission. Red contours are —4, —3 
(dashed), 3, 4, 5, and 6 (solid) times the rms of the map. The 
synthesized beam is shown in the bottom-right corner, and is 
1.74 X 1.64 arcsec 2 , with PA=— 89°. In both panels, the greyscale 
and black contours are the same as in Fig. [2}bottom. 
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Figure 4. Left column: spectra averaged in a region of ±20 arcsec 
with respect to the position of J041757. Right column: idem for 
J041757-NE (offset of (37,30) arcsec with respect to J041757). 
Note the excess of blueshifted emission at velocities in the range 
2-4 km s" 1 for CO (1-0). For the case of N 2 H+, J041757 was 
fitted with one Gaussian, while for J041757-NE we did fit the 
hypcrfine structure of N2H+. 



6.0-7.5 km s -1 range. We made line area maps by integrat- 
ing the CO emission in ranges of 1 km s _1 wide, from 4.5 to 
8.5 km s _1 , and found that for this velocity range CO traces 
the large-scale structure of the B211 filament (located to the 
south-west of J041757, see Fig.[T}, with strong intensity gra- 
dients and no clear features associated with J041757 (Fig. [5} 
top). This was expected because CO is optically thick and 
spread out over all the Taurus region (e. g., Goldsmith et al. 
2008; Narayanan et al. 2008). However, we fitted with gaus- 
sians the three velocity components seen in the CO spectrum 
and found, for the emission line at 5 km s _1 , an excess in the 
blueshifted emission with respect to the gaussian fttQ The 
line area map integrated in the range 2-4 km s _1 , where 
the blueshifted excess is seen, indicates that this emission is 
strongest around J041757, suggesting a possible relationship 
(Fig. Stop). 

Concerning 13 CO(1-0), well suited to trace the gas col- 
umn density in the region because is optically thinner than 
12 CO (1-0), the averaged spectra also reveal at least 3 veloc- 
ity components for both J041757 and J041757-NE. In Fig.[5} 
bottom, we show the line area maps of the 13 CO emission 
integrated between 4.5 and 8.5 km s _1 in steps of 1 km s _1 . 
Among the four mapped velocity intervals, the emission inte- 
grated between 6.5 and 7.5 km s _1 is more prominent at the 
position of J041757-NE, and has an extension to the south 
and another extension or stream to the southwest, passing 
through J041757. This could be a hint of the molecular gas 
associated with J041757-NE and J041757 coming from the 
same physical region. 

Finally, the N 2 H + (l-0) averaged spectra for J041757 
and J041757-NE reveal emission associated with J041757- 



8 An excess with respect to a Gaussian fit was also found with 
the Kitt Peak Telescope in 12 CO(1-0) between 2 and 4 km s -1 
(O. Morata, priv. commun.), with the integrated emission also 
concentrated on J041757. 



NE (see Fig|4]-bottom). We fitted the hyperfine structure of 
the N 2 H+(l-0) transition for J041757-NE using the pack- 
age CLASS of the GILD AS software, and obtained the out- 
put parameter At = 0.059 ± 0.006 K0, VLsr = 7.02 ± 
0.02 km s~\ Av = 0.42 ± 0.17 km s~\ and an opacity 
of the FlF'-FiF =23-12 transition of - 0.1, which corre- 
sponds to an excitation temperature of 3.4 K and a N2H 4 " 
column density of 2.5 x 10 11 cm -2 . This corresponds to a to- 
tal H2 column density of ~ 9 x 10 20 cm -2 (assuming a N2rl + 
abundance of (2-5) x 10" 10 , Pirogov et al. 2003, Tafalla et 
al. 2004), and a mass of ~ 3 Mj up (for a size of ~ 20 arcsec 
in diameter), similar to the lower limit of the mass derived 
from the submillimetre continuum emission (Table [4]). To- 
wards J041757, the strongest hyperfine component of the 
N2H + (l-0) transition is detected (at 5a), and a gaussian 
fit yields 7.04 ± 0.08 km s -1 , very similar to the velocity 
of the J041757-NE source. This suggests that both objects 
are associated with dense material (as the critical density of 
N2H + (l-0) is ~ 10 5 cm -3 ) possibly coming from the same 
region of Taurus, as indicated also by the morphology of 
the 13 CO emission integrated between 6.5 and 7.5 km s 
(Fig. [3-bottom). 



4 SPECTRAL ENERGY DISTRIBUTIONS 
(SEDS) 

Both J042118 and J041757 are Spitzer/IRAC sources emit- 
ting at 350 /im. In order to further constrain its na- 
ture, we checked the NASA Extragalactic Database (NED) , 
the National VLA Sky Survey (NVSS) at 21 cm, and 
JCMT/SCUBA catalog at 450 and 850 fim (di Francesco 
et al. 2008), and found no data for SCUBA, no significant 
emission in NVSS above 1 mjy beam -1 , and no registered 
object in NED. To look for further archival data we used 
the new capabilities of the Virtual Observatory SED Anal- 
izer (VOSA, Bayo et al. 2008; Bayo et al. 2012 in prep.). 
The tool queries virtual observatory compliant catalogs in 
the ultraviolet (GALEX, IUE HPDP photometry), optical 
(Tycho-2, CMC-14, Hauck et al. 1997, Mermilliod & Mer- 
milliod 1994, Sloan DSS release 7, and IPHAS), near, mid 
and far-infrared regimes (2MASS, IRAS, MSX6C, AKARI 
IRC and FIS, GLIMPSE, UKIDSS DR7, WISE preliminary 
release), and allows to fit several grids of models with ex- 
tinction as a free parameter. Appart from the data already 
presented in this work and in Barrado et al. (2009), we find 
UKIDSS counterparts for both sources but with unreliable 
photometry in the H band, and WISE counterparts at less 
than 1 arcsec from the Omega2000 position of J041757 (NIR 
component B) and the Spitzer/IRAC position of J042118. 
Since Sloan DSS did not detect any source at the position 
of J041757 and J042118, we give only upper limits at the 
corresponding frequencies. In the Appendix we list the flux 

9 In the CLASS package of the GILDAS software, the first output 
parameter of a fit to the hyperfine structure of a given transition 
is At = / (Jjj(T ex ) — Ji/(Tb g ))T, where / is the filling factor, 
assumed to be 1, J„(T) = hv/k/(e hv / kT - 1), with k being the 
Boltzmann constant, T the temperature and v the frequency; T ex 
is the excitation temperature, Tj, g the background temperature, 
and t is the opacity of the main hyperfine component, F^F'— 
F\F =23-12 for the case of N 2 H+(l-0). 
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Figure 5. Top: 12 CO(1-0) channel maps in the J041757 region. Rms is 0.16 K (or 0.23 K for the panel showing the emission integrated 
from 2-4 km s" 1 ). Bottom: 13 CO(1-0) channel maps in the J041757 region. Rms is 0.05 K. For all panels, colorscale is in main beam 
brightness temperature (K), the beam is ~ 22 arcsec, and stars mark the position of J041757 and J041757-NE. 



density measurements and upper limits (Table |A"T1 and 1X2)) . 
and in Fig. |6] we present the SEDs for both J042118 and 
J041757. 



5 DISCUSSION 

5.1 On the nature of J042118 

J042118 is undetected in Sloan DSS (data release 8) 
but detected in UKIDSS (data release 7, Fig. HJ, WISE, 
Spitzer/IRAC+MIPS and the CSO (see Table \Kty . In order 
to test if the J042118 SED (§ 0J is easily reproducible by 
a very extincted photospheric model (as expected if it was 
a background object seen through the molecular cloud), we 
used VOSA tool (§HJ) and found that even allowing extinc- 
tion values up to A v ~ 15 mag, it is not possible to get a 
reasonable fit (reduced \ 2 values of the order of hundreds) 
neither with Kurucz (Castelli et al 1997) nor with Lyon mod- 
els (Allard et al 1997; Baraffe et al 1997, 1998; Hauschildt 
et al 1999; Chabrier et al. 2000; Allard et al 2001; Allard et 
al 2003, 2007, 2009), suggesting that the object belongs to 
Taurus. 

With the aim of shedding light on the nature of J042118, 
we compared the CSO emission to the Herschel Archive im- 
ages at 250 and 160 »rrPl. which are the wavelengths for 
which the Herschel beam (17 arcsec and 11 arcsec) is more 
similar to the CSO beam (10 arcsec). While at 160 /im 



10 Herschel Archive data from the 'Herschel Gould Belt survey' 
Key Program (see Andre et al. 2010) were downloaded (observing 
identification numbers 1342202090 and 1342202254). The Her- 
schel/PACS(SPIRE) data at 160(250) fim were processed with 
HIPE(madMap) version 8 and Scanamorphos version 14 (Roussel 
2012). The Herschel data towards B211/L1495 will be published 
in a forthcoming paper by Palmeirim et al. (in prep.) as part of 
the 'Herschel Gould Belt survey' project. 
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Figure 6. Spectral Energy Distribution of J042118 (black line) 
and J041757, NIR component B (brown line), after Barrado et 
al. (2009). 



there is no significant emission towards J042118, the Her- 
schel/SPIRE 250 /im emission shows an extended large-scale 
structure elongated in the southwest-northeast direction (see 
Fig-O, with a subcondensation matching well the position of 
J042118, within the telescope positional uncertainties. The 
location of J042118 at the tip of the elongated structure, 
which seems to be part of the extended structures of the 
Taurus complex, favors its membership to Taurus. In this 
case, the submillimetre source detected with the CSO would 
be tracing a disc and/or small (< 5 arcsec or 700 AU of 
radius) envelope with a mass of ~ 1 Mj up (Table |4)|. The lu- 
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Figure 7. Submillimetre emission towards J042118. Greyscalc 
and black contours correspond to the Herschel/SPIRE 250 fim 
emission. Contours range from 0.06 to 0.12, increasing in steps of 
0.02 Jy beam" 1 (la is approximately 0.009 Jy beam -1 ). Herschel 
beam at 250 fim, shown in the bottom-left corner, is ~ 17 arcsec. 
White contours are the CSO 350 fim emission (same contours as 
in Fig. [2}, and CSO beam is ~ 10 arcsec. 



minosity (upper limit) for J042118 estimated from the SED 
(Table \M\ would be < 0.0023 Lq, and the SED shape is 
reminiscent of Class I(/II) young stellar objects, with a flat 
slope between the NIR and the submillimetre range. This 
can be quantified by estimating the bolometric temperature, 
of ~ 140 kF 1 !. typical of Class I young stellar objects (e.g., 
Chen et al. 1995). The fact that the submillimetre emission, 
as recovered by the CSO, is unresolved could be indicative 
of most of the envelope being accreted/dispersed, with the 
submillimetre emission arising only from a compact disc. 

5.2 On the nature of J041757 

The nature of J041757 is discussed in Barrado et al. (2009), 
who present series of color-color and color-magnitude dia- 
grams from Spitzer and NIR data to conclude that the most 
likely scenario for J041757 (for either component A and/or 
B) is that the object is galactic. Subsequent infrared spectro- 
scopic observations by Luhman & Mamajek (2010) indicate 
that component A is most likely a M2 background star, while 
no conclusion could be drawn for component B because no 
spectrum was shown for this component in that work. This, 
together with the fact that source B has a raising slope in the 
NIR SED matching well the FIR and submillimetre emission 
from Spitzer/MIPS and CSO (Barrado et al. 2009), indicates 
that the most likely counterpart of the submillimetre source 
in J041757 is the NIR source B, which in turn is associated 
with centimetre emission (§ I3.2.2[) . Hereafter we will refer to 
source B (identified in the NIR and in the centimetre range) 
as J041757-B. 

11 This number must be regarded with caution because the SED 
for this source is better sampled in the IR than in the submil- 
limetre and radio range, and because we treated all flux densities 
which are upper limits as detections (c.f. Kauffmann et al. 2008). 



Similarly to the case of J042118, we compared the 
CSO emission of J041757 to the Herschel Archive images at 
250 and 160 fim. See Fig. [5]to visualize such a comparison. 
The Herschel/SPIRE emission at 250 fim (black contours 
in Fig. |8]-top) is strongest at the position of J041757-NE, 
matching well the 350 fim emission observed with the CSO 
(white contours). From J041757-NE, there are two streams 
of emission, extending for about ~ 1 arcmin, one towards 
the south and the other towards the southwest and passing 
through J041757-B. The 250 fim emission towards the south- 
western side of the map is part of the B211 filament (Fig.[l|. 
We note the remarkable similarity between the emission at 
250 fim and the 13 CO(1-0) emission integrated between 6.5 
and 7.5 km s -1 (colorscale in Fig. |8]-top) . Compared to Her- 
schel data, the CSO is recovering the most compact (as 
we used the 'compact' option in the reduction process, see 
§ 12. 3[) and fainter emission (la noise of Herschel/SPIRE at 
250 fim: ~ 14 mjy beam -1 ; la noise of CSO at 350 fim: 
~ 6 mjy beam -1 ). Regarding the emission seen with Her- 
schel/PACS at 160 fim, this is again strongest at J041757- 
NE, and the northeast-southwestern stream splits up into 
one slightly resolved source coinciding, within the position 
uncertainties of the telescopes, with the 350 fim source at 
the position of J041757-B. This suggests that J041757-B is 
associated with the hottest dust of the stream linking the 
J041757-NE source and the B211 filament. Flux densities 
measured by Herschel at 160 and 250 fim are consistent with 
the flux densities given here as measured with the CSO, but 
we refrain from presenting a detailed study of the Herschel 
data in this work as this will be the subject of a forthcoming 
paper. 

5.2.1 J041 757-B: an extragalactic object? 

We consider first the possibility that J041757-B is an extra- 
galactic object, which would be aligned by chance with the 
submillimetre stream linking the J041757-NE source and the 
B211 filament. In Fig. [9] we show a comparison of the SED 
for J041757-B to the recent compilation of average SEDs 
of radio-loud and radio-quiet AGNs of Shang et al. (2011, 
which are in good agreement with the classical SEDs of Elvis 
et al. 1994). The SEDs are normalized at 12 fim (a good in- 
trinsic indicator for AGNs, e.g., Gandhi et al. 2009) to the 
flux of J041757-B, and the error bars are the standard devi- 
ation of the fluxes for Shang et al. sample, being represen- 
tative of the scatter in the AGN fluxes at each wavelength. 
From the figure, it is seen that J041757-B SED is not com- 
patible with the radio-loud AGN SEDs, mainly because of 
its important drop from MIR to radio wavelength! 1 ^!. In the 
visible range, J041757-B SED is much more obscured than 
AGN SEDs. In order to simulate an AGN extincted SED we 
reddened the AGN average SED for radio-quiet using the 
extinction law by Fitzpatrick (1999; improved by Indebe- 
touw et al (2005) in the infrared, and R v ~ 3.1), and with 



We note however that if we estimate the radio-loudness from 
the classical definition, Seem /^^ 00 ^ (following Kellermann et al. 
1989), we obtain ~ 700, which would classify J041757-B as a 
radio-loud. As explained in the text, this would be still compat- 
ible with a radio-quiet AGN if the extinction for J041757-B is 
important. 
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Figure 8. J041757: in both panels, colorscale is the 13 CO(l- 
0) area integrated in the range 6.5-7.5 km s" 1 (as in Fig. [5} and 
white contours are the CSO 350 fim emission (same contours as in 
Fig-O- Top: Black contours correspond to the Herschel/SPIRE 
250 (im emission. Contours range from 0.30 to 0.54, increasing in 
steps of 0.03 Jy beam -1 (lc is approximately 0.014 Jy beam" 1 ). 
Hcrschel beam at 250 um, shown in the bottom-right corner, is 
~ 17 arcsec ( 13 CO beam is ~ 22 arcsec). Bottom: Black con- 
tours correspond to the Herschel/PACS 160 (im emission. Con- 
tours range from 0.14 to 0.24, in steps of 0.014 Jy beam" 1 (la 
is approximately 0.018 Jy beam" 1 ). Herschel beam at 160 fim, 
shown in the bottom-right corner, is ~ 11 arcsec, similar to the 
CSO beam. 



A v ranging from 1 to 10 (see figure, thin blue lines). The re- 
sult is that the J041757 SED would be compatible with an 
AGN SED extincted by A v ~ 5 in the visible-IR range. How- 
ever, the J041757-B SED is not fully consistent with radio- 
quiet AGNs either, because of its important contribution of 
submillimetre (and radio) emission compared to MIR. This 
could be indicating that there is an excess of submillimetre 
emission in the direction of the galaxy, which could be ex- 
plained by the Taurus stream joining the J041757-NE source 
and the B211 filament. Nevertheless, there would be an ex- 
cess of submillimetre emission from the galaxy candidate it- 
self, as revealed by the image of Herschel/PACS at 160 fim, 
which shows that in the stream there is only emission at 
160 fim at the position of J041757-B (similar case for the 
CSO 350 fim emission, Fig. [8}. In addition, flat spectral in- 




log X (u.m) 

Figure 9. Spectral Energy Distribution of J041757-B (black) 
compared to average SEDs for radio-loud (thick red) and radio- 
quiet (thick blue) AGNs from Shang et al. (2011). Thin blue lines 
are the extincted radio-quiet SEDs with A v of 5 and 10 mag. 
The spectrum of J041757-B in the visible will be presented in a 
forthcoming paper (N. Huelamo et al., in prep.). 



dices in the centimetre range in AGNs are not usual, as cen- 
timetre emission from AGNs is typically produced by syn- 
chrotron emission from relativistic electrons accelerated in 
jets (e.g., Giroletti & Panessa 2009; Ibar et al. 2009; Shang 
et al. 2011; Randall et al. 2012). Furthermore, a preliminar 
spectrum of J041757-B in the visible range taken with Osiris 
on the GTC (Huelamo et al., in preparation) shows no broad 
emission lines, while the extincted SED in Fig. [9] indicates 
that even with A v > 5 still strong and broad lines are ex- 
pected for an AGN. If J041757-B was a galaxy it should 
be detected in forbidden emission lines in the MIR range 
([OIV], [NeV]: e. Goulding et al. 2009) and in hardest 
X-rays (> 10 keVl 13 !. which could be the final observational 
tests to definitively discard the nature of J041757-B as a 
background galaxy. 

5.2.2 J041 757-B: a background object behind Taurus? 

Similarly to J041757-A (Barrado et al. 2009, Luhman & Ma- 
majek 2010), J041757-B could be a background object be- 
hind the Taurus complex. However, if J041757-A is a back- 
ground star, the IR source mainly contributing to the FIR 
and submillimetre emission must be J041757-B. Therefore, 
whatever the true nature of J041757-B is, the object is em- 
bedded in considerable amounts of dust, and hence should 
be associated with molecular gas as well. This molecular gas, 
if J041757-B does not belong to Taurus, should appear as a 
velocity component in the line of sight different from Taurus 

13 In X-rays, we could only set an upper limit for J041757-B, 
which is not detected in the ROSAT All-Sky Survey (RASS, Voges 
et al. 1999, covering the range 0.1-2.4 keV with a typical limiting 
sensitivity of 10 -13 ergcm _2 s _1 , Anderson et al. 2003). 
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velocities, while in a 9000 km/s band we did not detect any 
other line apart from those of Taurus. If the object is asso- 
ciated with dust and gas, it should be a young stellar object 
(most likely in the Class II/III phase) or an evolved star. 
However, the centimetre emission for Class II/III objects 
is typically due to non-thermal gyrosynchrotron emission, 
which is not resolved (e.g., Andre et al. 1992, Forbrich et 
al. 2007) and is polarized, while J041757-B has a flat (ther- 
mal) spectral index, is slightly elongated in the southwest- 
northeast direction, and shows no polarized emission (we 
checked this in the dataset of 2008 observations and had no 
detection of polarized emission). In addition, the lack of X- 
ray emission (no ROSAT detection) is also uncommon for 
Class II/III objects (e.g., Giidel et al. 2007). 

The possibility that J041757 is an evolved stellar object, 
either a carbon star or a giant, is also considered. Luminous 
carbon stars are evolved stars with more carbon than oxygen 
in their atmosphere, characterized by red near-infrared col- 
ors and spectra with no H2O (Liebert et al. 2000). The JHK 
colors for J041757-B match well within the JHK color-color 
diagrams for these stars (e.g., Quanz et al. 2010). How- 
ever, we used the K s magnitude of J041757-B (Barrado et 
al. 2009) and Table 8 from Quanz et al. (2010) to estimate 
the distance and find > 160 kpc, larger than the currently 
known extent of our Galaxy. Since the emission in J041757 
is resolved in the submillimetre range, such a large distance 
would imply an unrealistic large size for the envelope of a 
carbon star of ~ 9 pc. In addition, the carbon star SEDs 
have the peak at around 12 /im, and longwards the flux de- 
creases strongly (e.g., Liebert et al. 2000), which is clearly 
not the case of J041757-B. Thus, the properties of J041757- 
B are not fully consistent with being either a Class II/III 
young stellar object, or an evolved star. In fact, by making 
use of the data compiled to build the SED (§ [4}, and us- 
ing the VOSA tool as outlined in § 15.11 we found that the 
J041757-B SED cannot be reproduced with a very extincted 
photospheric model, providing further support against the 
galactic background object scenario for J041757-B. 

5.2.3 J041 757-B: a proto-brown dwarf candidate in 
Taurus ? 

Finally, J041757-B could be an embedded object belonging 
to Taurus. This is strongly suggested by the discovery of a 
narrow (~ 16 arcsec width) stream joining the J041757-NE 
source and the B211 filament and passing exactly through 
the position of J041757-B. Since both J041757-NE and the 
stream are detected at around 7 km s _1 (in N2H + and/or 
13 CO), both structures belong to Taurus. 

If J041757-B belongs to Taurus, its low luminosity, of 
< 0.005 Z/q (using the values in Table IA2[I . would imply 
a proto-BD nature (Barrado et al. 2009). This would nat- 
urally explain the centimetre emission with a flat spectral 
index, since flat/positive spectral indices, indicative of ther- 
mal free- free emission, are typically observed in young stellar 
objects due to shocks in the jet (e.g., Anglada et al. 1998; 
Beltran et al. 2001; Avila, Rodriguez, & Curiel 2001; Pech 
et al. 2010; Scaife et al. 2011; Ward-Thompson et al. 2011), 
and have also been found for Very Low Luminosity Objects 
(e.g., Shirley et al. 2007, Scaife et al. 2011). In these cases, 
the gas is ionized in shocks produced by the jet driven by 
the nascent star, and the centimetre emission is usually elon- 



gated along the outflow axis. This would also easily explain 
the elongation seen at both 6 and 3.6 cm for J041757-B. And 
interestingly, the 12 CO (1-0) excess detected between 2 and 
4 km s is reminiscent of wings observed in outflows driven 
by low-mass young stellar objects. In addition, if J041757- 
B is a proto-BD, its bolometric luminosity and centimetre 
luminosity match well the relation found for young stellar 
objects driving radiojets, proposed first by Anglada (1995). 
Following the recent correlations of Shirley et al. (2007) 
and Scaife et al. (2011) for the low-mass case, the expected 
centimetre luminosity, given the bolometric luminosity of 
J041757-B, is 0.0002-0.007 mjykpc 2 , consistent with the 
measured centimetre luminosity of 0.004 ± 0.002 mjykpc 2 . 
Last, if J041757-B belongs to Taurus, the properties of the 
extended envelope where the NIR source is embedded, of 
about ~ 1000 AU of radius and about ~ 5 Mj up , coincide 
with the extrapolation of the trend presented by Kauffmann 
et al. (2008) in a radius versus mass diagram for small c2d 
Spitzer cores (see Fig. 4a of Kauffmann et al. 2008), and 
more recently by Sadavoy et al. (2010) and Schmalzl et al. 
(2010) for the Taurus region. 

Thus, if J041757-B belongs to Taurus, its properties 
can naturally be explained simply as a scaled-down version 
of low-mass young stellar objects. 

5.3 On the nature of J041757-NE: a 
pre-substellar core? 

J041757-NE is the strongest submillimetre source detected 
in the field of J041757 with the CSO, and has been confirmed 
by APEX, IRAM30m, and Herschel images. This source is 
clearly associated with Taurus, as shown by the N2H + and 
13 CO emission detected at 7 km s _1 (Figs. [4] and [5|. Since 
we estimated a total mass (compact -1-extended emission) for 
J041757-NE of ~ 100 Mj up (§E3), and no point IR source 
(up to 24 fj,m) has been found associated with the submil- 
limetre peak, J041757-NE could be harboring a true pre- 
substellar core, i.e., a starless core which will form a BD 
in the future. A first requisite for this to happen is that the 
submillimetre source should be gravitationally bound. Given 
the linewidth measured from the N2H + (l-0) hyperfine fit 
(§ 13.2, 3|l . and an average radius of 10 arcsec (Table|4j|, we es- 
timate a gravitational mass of 140 Mj up (following equation 
4 of Pound & Blitz 1993), similar to the total mass mea- 
sured for J041757-NE, suggesting that J041757-NE might 
be indeed a pre-substellar core. Interestingly, its total mass 
and radius follow the correlation found between mass and 
radius for protostellar/pre-stellar cores (e.g., Kauffmann et 
al. 2008, Sadavoy et al. 2010; Schmalzl et al. 2010). Also, 
the fact that this pre-substellar core candidate is found in 
B211, the portion of L1495 with smaller number of young 
stellar objects (e.g., Schmalzl et al. 2010) supports the idea 
suggested by Schmalzl et al. (2010) that this region is at the 
very beginning of the star (or BD) formation process. 

The number of pre-substellar core candidates in the 
literature is very scarce. A pioneer study searching for 
pre-substellar cores and proto-BDs was carried out by 
Pound & Blitz (1993, 1995), who found about 4 cores 
with substellar masses in the Taurus and Ophiucus clouds. 
However, the cores were at the detection limit and all but 
one were not gravitationally bound. In addition, the authors 
indicate that the candidate objects should be re-observed 
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deeply to rule out the possibility that the cores extend 
considerably, making them massive enough to overcome the 
substellar domain. More recent studies in Gould Belt clouds 
report samples of starless and pre-stellar cores, with some of 
them being substellar, but lack the kinematical information 
required to assess if they are gravitationally bound (e.g., 
Ward- Thompson et al. 2007; Davis et al. 2010; Sadavoy et 
al. 2010). Thus, J041757-NE is an excellent pre-substellar 
core candidate, and high angular resolution millimetre 
observations would be extremely useful to study possible 
infall motions within the core, which would definitely 
indicate that the object is on the way to form a BD. 



measured size and mass, we conclude that J041757-NE is 
likely gravitationally bound and thus is a good pre-substellar 
core candidate. 

Overall, this work shows the power of submillimetre in- 
struments to trace extremely low masses and thus to reveal 
substellar objects. We found evidence of two proto-BD can- 
didates and one pre-substellar core candidate which have 
scaled-down properties of low- mass young stellar objects and 
which are embedded in large-scale structures, all this favor- 
ing the in-situ scenario for BD formation in this portion of 
the Taurus complex. 



6 CONCLUSIONS 

We present the results of a multiwavelength search for 
pre-substellar cores and proto-BD candidates in the B213- 
L1495 clouds of Taurus. The studied sample of 12 candi- 
dates, selected from Spitzer/IRAC photometry of the re- 
gion (Barrado et al. 2009), was observed at 1.2 mm with the 
IRAM 30m/MAMBO-II bolometer. The two most promising 
objects, J042118 and J041757, were subsequently observed 
at 350 /im with the CSO. Additional observations were per- 
formed toward J041757, to further constrain its nature. Our 
main conclusions can be summarized as follows: 

(i) None of the 12 Spitzer/IRAC selected sources observed 
with PRAM 30m/MAMBO-II were detected at a level > 4a, 
and two sources, J042118 and J041757 had fluxes at the 
2a level. We set upper limits to the envelope masses of < 
3 Mjup. 

(ii) CSO 350 /im observations towards J042118 reveal one 
compact submillimetre source associated with the selected 
Spitzer/IRAC source, with a size of < 10 arcsec ~ 1400 AU 
(at the Taurus distance), and a mass of ~ 1 Mj up . A prelim- 
inar analysis of the Spectral Energy Distribution indicates a 
luminosity of about 0.002 Lq, and a bolometric temperature 
around 140 K, typical of Class I/II young stellar objects. 

(iii) CSO 350 /im observations towards J041757 reveal one 
extended source associated with the selected Spitzer/IRAC 
source, of about ~ 16 arcsec or ~ 2000 AU of size, and 
a mass of ~ 5 Mj up . Additionally, we detect one sligthly 
resolved source at 3.6 and 6 cm exactly at the position of 
the infrared counterpart 'B' of J041757, called J041757-B, 
with a fiat spectral index. J041757-B is found at the centre 
of a stream of dust and gas, extending from the B211 main 
filament to the northeast, mainly detected with 13 CO(1-0) 
at a velocity of ~ 7 km s _1 . N2H + (l-0) emission is detected 
at the same velocity at 5a towards J041757-B. All this seems 
to indicate that J041757-B is probably a proto-BD in the 
Class 0/1 stage associated with a radiojet and an extended 
envelope, whose properties seem to be a scaled-down version 
of low-mass protostars. 

(iv) We discovered a partially extended submillimetre 
source about 40 arcsec to the northeast of J041757-B, of 
20 arcsec ~ 3000 AU of size and total mass of ~ 100 Mj up , 
called J041757-NE, and with no centimetre emission neither 
optical/infrared (up to 24 /im) emission associated. This in- 
dicates that J041757-NE seems to be starless. J041757-NE 
is well detected in 13 CO and N2H + , and is found also at 
~ 7 km s _1 . With the measured linewidth of N2H + , and its 
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APPENDIX A: PHOTOMETRY FOR J042118 
AND J041757-B 

In this Appendix we give the photometry for J042118 and 
J041757-B, from visible wavelengths up to the upper limit 
set at 1.2 mm with MAMBO-II. To compile the photome- 
try we used the Virtual Observatory Spectral Energy Dis- 
tribution Analyser (VOSA, Bayo et al. 2008). See Hewett 
et al. (2006) and Lucas et al. (2008) for a review of the 
United Kingdom Infrared Digital Sky Survey. For the case 
of J041757-B, we measured again the flux at 24 /xm to take 
into account extended emission (neglected in our previous 
measurement given in Barrado et al. 2009). 



Table Al. Photometry for J042118 
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Table A2. Photometry for J041757-B 
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